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GENERALIZED SOLUTIONS 
TO LINEARIZED EQUATIONS 
OF THERMOELASTIC SOLID 
AND VISCOUS THERMOFLUID 

ANVARBEK M. MEIRMANOV, SERGEY A. SAZHENKOV 



Abstract. Within the framework of continuum mechanics, the full descrip- 
tion of joint motion of elastic bodies and compressible viscous fluids with tak- 
ing into account thermal effects is given by the system consisting of the mass, 
momentum, and energy balance equations, the first and the second laws of 
thermodynamics, and an additional set of thermomechanical state laws. The 
present paper is devoted to the investigation of this system. Assuming that 
variations of the physical characteristics of the thermomechanical system of 
the fluid and the solid are small about some rest state, we derive the linearized 
non-stationary dynamical model, prove its well-posedness, establish additional 
refined global integral bounds for solutions, and further deduce the linearized 
incompressible models and models incorporating absolutely rigid skeleton, as 
asymptotic limits. 



1. Introduction 

We are interested in proposing of a mathematical description of small pertur- 
bations in the thermomechanical system consisting of interacting elastic solid and 
viscous compressible fluid. In our study we aim to integrate purely mechanical, 
thermodynamical, and heat transfer effects altogether under one umbrella. Within 
such unified approach, a thermoconductive elastic body may be named a thermoe- 
lastic solid. Also we use the term thermo fluids, which has been introduced quite 
recently for a subject that analyzes systems and processes involved in energy, vari- 
ous forms of energy, and transfer of energy in fluids [111 Sec. 1.1]. 

The basic mathematical concept of reciprocal motion of thermoelastic solids 
and viscous thermofluids incorporates the classical conservation laws of continuum 
mechanics, the first and the second laws of thermodynamics, and a set of state laws 
specifying individual thermomechanical behavior of the components of media. In 
this article this concept is called Model 0\ and is stated in the beginning of Sec. 
[2] It is quite universal and spans a large variety of different phenomena in nature 
and technology. (A rather general relevant observation may be found, for example, 
in [TTJ Sees. 1. 5-1.7, VI].) At the same time this model is very complex and highly 
nonlinear. Therefore, some physically reasonable simplifications are necessary in 
view of further applications to natural problems and in engineering. 
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In investigation of small perturbations it is suitable to simplify Model 0\ by 
implementing the classical formalism of linearization about a rest state [3j Sec. 
V.7]. As the result, the linearized model arises, whose core consists of the heat 
equation coupled with the non-stationary compressible Stokes system in the fluid 
phase and the system of the wave equations (in the elasticity theory also called 
Lame's equations) in the solid phase. The linearization procedure and the precise 
formulations of the resulting Model O2 and its dimensionless version Model A are 
outlined in Sees. and [3J 

Also in Sec. [3Jwe introduce the notion of generalized solutions of initial-boundary 
value problem for Model A and formulate the first main result of this article - 
Theorem 13.41 on existence and uniqueness of solutions to Model A. The proof of 
this theorem relies on classical methods in the theory of generalized solutions of 
equations of mathematical physics and is fulfilled in Sees. H]and[5l 

After this, we are interested in studying of the limiting regimes in Model A, 
arising as some coefficients grow infinitely. To this end, in Sees. [6] and [7] additional 
global integral bounds are established for the pressure distributions (see Theorem 
16. 3p and for the deformation tensor in the solid phase (see Theorem I7.3[) . These 
bounds constitute the second main result of this article. With their help, in Sees. 
[5] and [5] Models Bl, B2, and B3 of incompressible media and Models CI and C2 of 
viscous thermofluid contained in an absolutely rigid heat-conducting skeleton are 
established as respective incompressibility and solidification asymptotic limits of 
Model A (see Theorems 18.11 l9Tl and l9.4j) . These models are the third main result 
of this article. 



2. Basic Nonlinear Formulation and Linearization 

Let Q be an open bounded set in K 3 with a smooth boundary. From now on we 
assume that f2 is the cube with the side of a size Lq, i.e., = (0, Lq) 3 . Next, assume 
that at time t = O is occupied by a solid component f2 s and by a fluid component 

flf = CI \ fi s such that the interface between the components T = dfl / fl dfl s is a 
rather smooth surface or a finite union of such surfaces. The both solid and fluid 
phases obey the fundamental conservation laws, which have the following forms in 
Lagrangian variables x and t [3J Sec. V.5], [HI Sec. 10]. The balance of mass 
equation is 

po = pJ, x in flf or £l s , t > 0, (2.1) 
the balance of momentum equation is 

= ^div x (TP) - JTPTTMiv^T*)- 1 + p Q F, xmn f ortt s , t > 0, (2.2) 
and the balance of energy equation is 

p ^ = -Jdiv^T" 1 ^ + Jq ■ div^TT 1 • ./, : : • 
at ot 

q = -x(T*)- 1 V a '!?, ccm%orO s , t > 0, (2.3) 

where 

T = V x f , r-x = w, J = dctT, 2E = T*T-I. (2.4) 
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The thermodynamical state of the both phases is governed by the first law of 
thermodynamics 

Us = dU + pd- (2.5) 
P 

and a thermodynamical state equation 

U-ds = F{p,$). (2.6) 

We postulate that thermomechanical behavior in the solid phase is described 
by the Duhamel-Neumann law of linear thermoelasticity, which is consistent with 
thermodynamical relations (|2.5[) and (|2.6p Sec. 1.6], [T^l Sec. 10]: 

P P s = f -Po - VJs(V -$*)+(v- I A div x w] I + 2XB(x, w), 



3 



x e Q s , t>0, (2.7) 



and that thermomechanical behavior in the fluid phase is described by the Stokes 
state equation, which has the following form in Lagrangian variables [3J Sec. V.5], 
P2 Sec. 6]: 

P t/ P/ = (_p + (, _ 2 r )tr ((T*)" 1 ^ T- 1 ))! + 2 / z((T*)- 1 ^T^ 

x GO/, t>0. (2.8) 

In Eqs. (|2.ip - (|2.8p w, p, p, U, and s are unknown displacement, temper- 
ature, density, pressure, specific intrinsic energy, and entropy, respectively. The 
medium under consideration is a two-parameter thermomechanical system. Fur- 
ther we refer $ and p to as the independent thermodynamical parameters. The 
triple r — (ri,r 2 ,r 3 ) is the set of spatial Eulerian coordinates of particles of solid 
or fluid. The vector q is the heat flux. The second equation in (|2.3[) is the Fourier 
law for the heat flux. 

T is the distortion tensor. T', T _1 , and (T') _1 are its conjugate, inverse and 
inverse conjugate, respectively. E is the deformation tensor. It is connected with 
the displacement vector by the formula 

2E = {V x wf + V x w + (V x wYV x w. 

B(x, to) is the symmetric part of the gradient W x w, i.e., 

2B(x, w) = {V x wY + V x w. 

P is the stress tensor. Note that, due to the identity [3l Sec. V.4], [T2l Sec. 6] 

5 = T*©(r, v)T, 
dt v ' ' ' 

where U>(r,v) = (l/2)((V,-7;) t + V r u) is the deformation rate tensor in Eulerian 
variables r and t, and v is the velocity of particles of fluid or solid, Eq. (|2.8p 
reduces in Eulerian variables to the well-known form 

P f = (-p+ (v-^div r v)I+2pM(r,v), r € t>0. 

Coefficients x, p, r), A, and 7 S are given. They are a thermal conductivity, 
bulk and shear viscosity coefficients of the fluid, bulk and shear elastic modules of 
the solid, and a thermal extension of the solid, respectively. In general, they may 
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depend on the thermodynamical parameters i9 and p. In line with the thermody- 
namics fundamentals, we have [El Chap. 5], [9l Chap. 1] 

2 2 

x, v, /i, rj, A, 7 S > 0, ^ > -n, t] > -A. (2.9) 

Constant positive coefficients and po are given temperature of some rest state 
and atmosphere pressure, respectively. 

Functions po, F, ^f, and T are given. They are an initial distribution of density, a 
density of mass distributed forces, a volumetric density of exterior heat application, 
and a specific free energy, respectively. Thermodynamical behavior in the solid 
and the fluid is different, in general. Therefore we are given, in fact, two distinct 
functions T s and Tf such that 

' Tf{p,&) if xen f , t> 0, 
F a \p,$) if xefl s , t>0 

and two distinct functions \& s and f such that 

*/(jc,t) if xeftf, t > 0, 
^ s (x,t) if x e Q, s , t>0. 

Similarly, we have two distinct coefficients of thermoconductivity, *c s in the solid 
phase and Hf in the liquid phase. 

In terms of the indicator function of the fluid phase 



*(a,t) 



we may write 





= %FfM + 


(I-x)^(p,^) Vxe!J, 


t > 0, 


(2.11) 




= X*/(x,*) + 




* > 0, 


(2.12) 


x(p,0) 


= X*f{p,tf) + 


(l-x)x,(p,0) Viefi, 


t > 0. 


(2.13) 


Analogously, 












P = XF/ + (1 


-x)P a Vccefi, t>0. 




(2.14) 



Interactions between the fluid and the solid are governed by the classical condi- 
tions on discontinuity surfaces [31 Sec. II. 3], [TJl Sec. 12], [T31 Sec. 2]. In order 
to state these conditions, we introduce the following notation. For any x £E T and 
for any function <p(x), continuous in the interior of O s and in the interior of Qf, 
denote 

f(s)(xo)= lim ip(x), tpm(x )= lim ip(x), 

X — » X|) x — * JBQ 

[<p](x ) = tp( s ){xo) - V(/)(^o)- (2.15) 

Remark 2.1. Clearly, T is immovable in Lagrangian coordinates in the sense that 
its any parametrization in Lagrangian variables does not depend on t. 

Thus, in terms of notation (|2.15|) we write down the conditions on T as follows: 
the continuity of temperature 

[tf] = 0, x e r, t > 0, (2.16) 

the continuity of displacement 

[w] = 0, x e r, t > 0, (2.17) 
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the continuity of normal stress 

[Fn] = 0, x o er,t>0, (2.18) 
and the continuity of the normal heat flux 

[*r(T t ) _1 V a; i? • n] = 0, x eT, t>0. (2.19) 

In Eqs. (|2.18p - (|2.19p n(x ) is the unit normal to T at a point x e T. We 
suppose that n is pointing into flf. 

The conditions of continuity of temperature and displacement on T reflect the 
local thermodynamical equilibrium and no-slip effect on the interface, respectively. 
Besides, the no-slip condition (|2.17[) includes the mass conservation law, applied to 
r, and manifests that T is a contact discontinuity, which means that the solid and 
the fluid do not exchange particles. The conditions of continuity of normal stress 
and normal heat flux are the respective momentum and energy conservation laws 

on r. 

Eqs. (|2.ip - (|2.8p and (|2.16[) (|2.19[) constitute the closed nonlinear model of joint 
motion of thermoelastic solid and viscous thermofluid. By the standard procedure 
[3j Sec 8.1], thermodynamical parameters s, p 7 and U may be expressed in terms 
of p and ■& from Eqs. (|2.5p and (|2.6|) by virtue of the free energy T: 

fl = -^G>,0), p = p 2 T' p { P J), U=-W#(p,ti)+F(p,ti). (2.20) 

Insert Eqs. (|2~20|) into Eqs. ([2~3jl and (|2~8| . and discard Eqs. (|2~5|) and (f2T6|) 

to get the equivalent system of equations and interface conditions, which we name 
Model 0i. 

Now simplify Model 0\, applying the classical formal procedure of linearization 
[31 Sec. V.7], [Ml Sec. 2.2-2.3] to it. Assume that p s and pf are mean constant 
densities of the solid and the liquid at rest. Suppose that the temperature in this 
rest state is equal to the constant temperature which was introduced in Eq. 
(pTfp . Denote 

c s =J r s(Ps,'&*), Csp = F' sp (Ps,fi*), c s # =^(p s ,i?*), 

Cspp = F'sppiPs, 1 ®*), C spi} = fJpviPs^*), CsOO = J r %(Ps,'&*), 

C*f=Ff(pf,&*), Cfp=J : ' fp (pf,&*), c /l9 =J : ' f #(p f ,$*), 

c f PP = FfppiPf^*), c fP® = FfpoiPf^*), c m = FfwiPf,^*). (2.21) 

Expand Eqs. (|2.20p for p and U in Taylor's series and skip terms of orders higher 
than one to get the linearized expressions 

p{p, d) = Cfpp) + {2c fp p f + c fpp p)){p - p f ) + Cfpvp)^ - i?*), (2.22) 

U s {p,'&) = -c slJ $* + c* + (c sp - c sp #??*)(p - p s ) - c si} #'3*('d - (2.23) 

U f (p, i?) = -Cftf, +c* f + (c fp - c /t #4 m ){p - p f ) - c /lM i?*(tf - (2.24) 

Remark 2.2. In Eqs. (|2.22p and (|2.24|) notice that p(pf, t?*) = po due to interface 
condition (|2.18p . Also notice that from the physical point of view it is observed 
for the most of continuous media that p(p, $) is increasing in p and $ and that 
Uf(p,'d) and U s (p,'d) are increasing in d and decreasing in p. Thus Cf p $ > 0, 
2c /p + c fppPf > °' c s-m < °' °sp - c-sptifi* < 0, Cf&g < 0, and c/ p - Cf p &d < 0. 

Substituting Eqs. (|2.22p - (|2.24p into Model Oi, expanding the terms of thus ob- 
tained equations in Taylor's series in p, w, and derivatives of w with respect to 
x and t about the rest state, and discarding the terms of orders higher than one, 
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we arrive at the following linearized dynamical model of the system thermoelastic 
body - compressible viscous thermofluid, which we call further Model GV 

Statement of Model 2 - The motion in the solid phase is governed by the 
classical equations of the linear thermoelasticity theory: 

Ps ~dW = dWxFs + PsF > xen s , t>0, (2.25a) 

d-d d 
-p 5 c sM i?*— = x s A x d - i)^ s r]—div x w + * s , x & fl s , t > 0, (2.25b) 

P,= (-Po - 7.f(0 - *0 + (v-^dW xW y + 2XD(x,w) 7 xen s ,t>0, 

(2.25c) 

p = p s (l - div x w), xett s ,t>0. (2.25d) 

The motion of the liquid phase is described by the linearized classical model of 
liquid and gases: 

d 2 w 

Pf^2- = d[v * ¥ f + Pf F > xeQ f ,t>0, (2.25e) 

Qft Q 

-pfCfwd*— = x f A x i9 - dvCfpvpj— div x w + */, xeflf, t > 0,(2.25f) 

¥ f = (-p - CfpdPfifl - ■#*) + (2c f p + c fpp pf)p 2 f div x w 

+ (i/-|/i)div I ^)l + 2/iD(a;,^) ) x e Sl f , t > 0, (2.25g) 

p = p f (l - div x w), xetlf, t>0. (2.25h) 

The linearized conditions on the interface between the solid and the fluid are 

[tf] = 0, x e T, t > 0, (2.25i) 

[to] =0, x„er,f>0, (2.25j) 

(P s ( 1 ? (s) , W(s) )-P / (% ) , W(/) ))n = 0, x Q eT,t>0, (2.25k) 

(x s V x tf (s) - x/V x %)) -n = 0, xo e T, t > 0. (2.251) 

Remark 2.3. Transforming Model 0\ to Euler variables (r,i) and applying the 
same linearization formalism as above, one arrives exactly at Model Oi- Therefore 
in linearized setting the Lagrange and Euler descriptions of the thermomechanical 
system under consideration coincide. 

3. Problem Formulation 
and Statement of the Existence & Uniqueness Theorem 

In line with later asymptotic analysis in the present and forthcoming papers, we 
bring Model O2 to a dimensionless form and absorb the interface conditions on T 
in the equations by introducing a uniform description of the both phases. 

More precisely, choose the diameter L of the domain fi, a characteristic du- 
ration of physical processes r , acceleration of gravity g, atmosphere pressure po, 
mean density of air po at the temperature 273 K under atmosphere pressure, and 
the temperature difference $0 between the thawing and freezing points of water 
under atmosphere pressure as characteristic scales of length, time, density of mass 
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distributed forces, pressure, density of matter, and temperature, respectively, and 
denote 9 = d — 

Next, introduce the dimensionless variables (with primes) by the formulas 
x = Lqx' , t = Tot' , w = Low' , 



P = PoP 



P = Pop 



(3.1) 



the dimensionless vector of distributed mass forces, volumetric densities of exterior 
heat application, thermal conductivity coefficients, and mean densities of the solid 
and the fluid at rest (all with primes), respectively, by the formulas 



gF', * s 



ro^o 



Po$* 

T~0$0 



_ Llpoi)* 



'0 To^l 
and the dimensionless ratios by the formulas 



Poih , 



<f, Ps = PoP s , Pf 



PoPp 



(3.2) 
(3.3) 



a,-, 



4 



n 



l 

Pa 
JoPo 
Po 

TOPo' 



Cpf 



c 

2A 

Po ' 

Po ' „ 
CjmPfdp 

Po 



1 / 2 
i»n v 3 



TOPO 

7oc 2 , 
—Pf 



c; 



ii 



010 f 



Po ' 

Cs&tfPs$o 
PO 



(3.4) 



Quantity Co is the speed of sound 



where c = yj2cf p pf + C} pp pj and 70 = 7/5. 
in air at the temperature 273 K under atmosphere pressure. Quantity c is the 
speed of sound in the considered fluid at the temperature under atmosphere 
pressure. Dimensionless constant 70 is the ratio of specific heats (in other terms, the 
polytropic exponent) for air at the temperature 273 K under atmosphere pressure 
[Til Appendix 4], [Jj Sec. 2.4]. 



Remark 3.1. On the strength of Eqs. (12. 9| and Remark 12.21 all dimensionless 
constants on the left-hand sides of relations (|3.4|) are positive. 



Now, first shift the pressure scale on the constant value po so that the stress 
tensors (|2.25cj) and ( |2.25gP become vanishing at the rest state. Next, multiply Eqs. 
(I2.25aj) and (|2~2l3e)) by L/p , Eqs. (|2.25b|) and ([2~25f]) by t ^ /(po^*): and divide 
Eqs. (|2.25c[) and Q2.25gP by po- After this, substitute expressions (|3. ![) — (|3.4|) into 
the resulting equations and then omit primes. Thus, Model O2 is brought to a 
dimensionless form. Finally, use the notation (I2.10p - (I2.14[) and additionally set 



P = XPf + C 1 - X)Ps, 
c P = XCpf + (1 - X)c ps , 



ae = xaef + (1 - x)ote s , 
x = x>ef + (1 - x)xs 



(3.5) 



and introduce the dimensionless pressures p, q, and ir in order to wrap the dimen- 
sionless Model O2 into the following form. 



Statement of Model A. In the space-time cylinder Q = il x (0,T), where 
il = (0, l) 3 and T = const > 0, it is necessary to find a displacement vector w, a 
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= x 



temperature distribution 9, and distributions of pressures p, q, and ir, which satisfy 
the equations 

_d 2 w 

a rP-^r = div^P + a F pF, (3.6a) 
ot 

( X, ^dl)) +( 1 -x)(-^ + a x ^(x,w)) -O001, (3.6b) 

39 d 
tv ~di = diVx (^ x °} ~ a g —div x w + (3.6c) 

P + X a pdiv x w = 0, (3.6d) 

q = p^ j-, (3.6e) 

a p ot 

n + (1 — X) a rid\v x w — 0. (3.6f) 

We endow Model A with initial data 

dw 

w\ t=0 = w , — \t=o = v , t=o = o , xefl (3.6g) 
ot 

and homogeneous boundary conditions 

w = 0, 9 = 0, x e dQ, t>0. (3.6h) 

Remark 3.2. On the strength of classical theory of conservation laws of continuum 
mechanics [3 Sec. II.3], Eqs. |3?6a| - (|3.6f|) yield interface conditions (|2.25i|l - (|2.251|) 
(in the dimensionless form). Thus Eqs. (|3.6ap — (]3.6fj) are equivalent to Model O^- 

Generalized solutions of Model A are understood in the following sense. 

Definition 3.3. Five functions (w,9,p,q,n) are called a generalized solution of 
Model A if they satisfy the regularity conditions 

dw . . m / dw 1 
-,B( X ,w), X B( X ,- 

boundary conditions (|3.6hp in the trace sense, equations (|3.6d|) — (|3.6f|) a.e. in Q, 
and the integral equalities 



to, B(x, w), x©(x, ^) , div x w, 9, V x 9 e L 2 (Q), (3.7) 



f ( _dw dip \ f 

J [a T p— ■ — -P :X7 x ip + a F pF ■ ip jdxdt + J a T pv ■ <p\ t=0 dx = (3.8) 

for all smooth ip = (p(x,t) such that ip\on — <p\t=T — and 

c p 9 + agdiv x w )ip\ t=0 dx = (3.9) 



for all smooth -ip = ^P(x,t) such that ip\dQ = ip\t=T = 0. 

The first main result of the paper is the following theorem on existence and 
uniqueness of solutions to Model A. 

o 

Theorem 3.4. Whenever w eW% (A), v , 9 € L 2 (Q), and F, * e L 2 (Q), Model 
A has a unique generalized solution (w, 9, p, q, 7r) in the sense o f Definition 1 3. 3\ 
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4. The Energy Estimate and Uniqueness of Solution of Model A 

Construction of the energy estimate is based on introducing the alternative equiv- 
alent definition of generalized solutions of Model A and on a special choice of test 
functions in the integral equalities in this definition. 

Namely, we state: 

Definition 4.1. Five functions (w,9,p,q,ir) are called a generalized solution of 
Model A if they satisfy the regularity conditions 

dw m , , / dw - 

w, 



dt 



, B(x, w), xp(x, , div x w, 9, V x 9 e L 2 (Q), (4.1) 



functions t 1— > —— (t), t i— > c p 9(t) + agdiv x w(t) are weakly 



dw 
~dt 

continuous mappings of the interval [0, T] into L 2 (f2), (4-2) 

boundary conditions (|3.6hp in the trace sense, equations (|3.6d[) — (|3.6fj) a.e. in Q, 
and the integral equalities 

J J ( a rP^f ■ ^ - p : V x v + a F pF ■ cp^dxdt 

f dw f 

= / a rP^rr( x , T ) ■ <p(x,T)dx - / a T pv ■ ip\ t = Q dx, Vre[0,T] (4.3) 
Jn ot Jn 

for all smooth <p = tp(x,t) such that <p\gn — and 

f [ (c„9 + c(0div x w] — kV x 9 ■ V x ijj + ^ip dxdt 
Jo 7o LV J dt J 

= / (c p 9(x, t) + a s div x w(x,T))if;(x,T)dx - / (c p 9 + agdiv x w )ijj\t = odx, 
Jn Jn 

Vre[0,T] (4.4) 

for all smooth ip = ip(x,t) such that ifi\dn = 0. 

Remark 4.2. Definitions 13.31 and 14.11 are equivalent. The fact that a generalized 
solution in the sense of Definition 14.11 is a generalized solution in the sense of 
Definition ^. 31 is quite obvious. The inverse proposition is true thanks to the simple 
standard considerations. Its justification can be fulfilled similarly to, for example, 
[D Sec. III.l]. 

Now let us follow the track of considerations of [TU1 Chap. 2, Sec. 5.2]. 

Fix arbitrary t*,t«* € (0,t), t* < r**. Take a continuous piece- wise linear 
function on [0, r] such that <j> m (t) = 1 if r* + (2/m) < t < t** — (2/m) and 
<t>m(t) = if t > r** — (1/m) and < < t* + (1/m). Take a regularizing sequence 
Lu n £ C£°(R) such that 

1 1 

n n 



uj n {t) = uj n (-t), w n (t) > 0, / uj n {t)dt = 1, suppw„ C 
For n > 2m set 

dw 
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where the asterisk* means the integral convolution in R, and substitute for <p and 
i/j into Eqs. (|4.3|) and (14. 4j) . respectively. Clearly, this choice of test functions is 
valid due to regularity properties (|4.ip and (|4.2[) . 

Insert Eqs. (j3.6ci[) — (|3.6f|) into Eq. (|4.3[) (with <p = <p„ m ) and then sum the result 
with Eq. (|4.4[) (with V- 1 = ipmn)- I n thus obtained equality represent 



o ./n 



a00div x ip mn dxdt 



o Jn 



a g (((f) m 9) * u n ) —div x ((<f> m w) *uj n )dxdt 



o Jn 



ae{{(l)mO) * u„)div x ^-^wj *uj^dxdt (4.6) 



and 



o 



agdrv x w — - — axclt 
dt 



o Jn 



agdiv x w((4> 



-dxdt 



o ./n 



+ / / a e div a; ((^ m w) ({<t> m Q) * u n )dxdt. (4.7) 



Integrate by parts with respect to t in the last summand of Eq. (|4. T|) and cancel 
similar terms. 

Applying the arguments of Chap. 2, Sec. 5.2], after some technical trans- 
formations and passage to the limit as n / oo, m / oo, t* \ 0, and t*„ /* t, 
successively, we finally arrive at the energy identity, as follows: 



1 

2°* 



Oil) 2 j 



2 

2, a ' 



-apHxdiv^T) 



2 

2, a 



dw 2 
Xd 1Vx — 



2 

dw\ I 2 



- 



dt ) 



2,f2x(0,r) 

ia T ||V^^o||2,o + TTQ^IK 1 - x) div xW>o||lfi + ^Mxdiv^ttfoIlL 



2,fix(0,r) 

1 

— ( 

2 



x°\\2,nx(o,T) 



1 



1 



-a A ||(l-xM^,^o)IHo + oll^ol 



2, SI 



o Jn 



{ctFpF 



dw 
~dt 



^6) dxdt, Vre[0,T]. 



Discarding all the terms on the left-hand side except for the first and the fifth 
ones and applying the Cauchy-Schwartz inequality on the right-hand side, we get 



r=dw . 2 

~dt 



1 



2.12 



2. SI 



< 7; a r\\Vpi>o\\in + ollV^p^oll^n 



T 1 



^dw , , 



dt 



T 1 



emlndt+^W^pFW 



l2,fJx(0,- 



2min{c p /, c ps } 



2.Q 



2,Ox(0,r) 



1 



+ ^IK 1 - x)div x w ||2,n + ^"pllxdiv^ol^ 



+ 2 a ^\0--Xp(x ) w )\\l n) Vre[0,T]. (4. 
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Applying Grownwall's lemma to this inequality, we conclude that 



1 



2,nx(0,r) 2 



O II \Ap^ll2,nx(0,r) 



< 



T . „2 



|^llvWll,nx(o,t) + 2wm{c pf ,c p .} m!l > a *« ) ' e > 



e T ~ t dt 



+ e 7 



rarllv^wolll.n + ollv^P^llIn 



+ (e r - 1) IK 1 ~~ X)^ v xw a \\l n + ~a p \\xdiv x w \\2 iQ 

+ ~ax\\{l-xM*>Vo)\\l a ], Vr€[0,T]. (4.9) 

Combining energy identity (|4.8|l with estimates (|4.8jl and (|4.9|l . we finally estab- 
lish the following. 

o 

Proposition 4.3. Let F, * £ L 2 (Q), w eVt^ 1 (Q), and v ,9 € L 2 (Q). Assume 
that five functions (w, 9,p, q, w) are a generalized solution of Model A. Then w and 
9 satisfy the energy estimate 



- ess sup a 7 
' te[o,r] 



Quo ^ I 



te[o,T] 



+ ^ ess sup a v ||(1 - x)<ii«a;«>(*)||| n + \ess supa p \\xdiv x w(t)\\ 2 2n 
z te[o,r] ^ te[o,T] 



1 



+ -esssupaA||(l-xM^^(*))ll2,n + °v 



te[o,r] 



9k; 2 
ydiv x -— 
A dt 



2,f2x(0,r) 



\1 



2,f2x(0,r) 



^\\VpF\\la xm + 2m i n { Cp/)Cp5 } ll * l| 2^x(o,t) 



e r + 1 



Hrllv^Wol 



2 

2. SI 



c p y O|| 2 ,f2 



a^lK 1 - x)div x w \\l n + apWxdiVxWoWl n 



+ a x \\(l- X P(x ) w )\\l Q ='c en (T), Vre[0,T]. (4.10) 



Corollary 4.4. The displacement w belongs to the space L°°(0, T;W2 (Q)) and 
admits the bound 



ess sup \\w(t)\\ w i m < C(C en (T),C k (n)), 
te[o,T] 



(4.11) 



where Cfc(fi) depends only on geometry of <9fi. 
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Proof. On the strength of the Newton-Leibnitz formula, 

dw 



x(x)H)(x, w (r)) : G(x)dx = / / X {x)B[x, — (t) : &(x)dxdt 

Jo Jn \ at J 

+ [ x(xp(x, w ) : G(x)dx VG e (I 2 (fi)) 3x3 , Vre[0,T]. (4.12) 



Substitute Gr(x) = H)(x, w(t)) and use inequality (|4.8j) and the Cauchy-Schwartz 
inequality with some positive £o and e± to get 



esssup / x|D(x, w(t))\\ 2 dx 
te[o,r] Jn 



< 



o Jn 



1 



X 



2e 

r l 

2Ti 



dw 



^-(t)) +f2|D ( ,-. tt ,(r)l 



dxdt 



Kx,w )\ 2 + ^\B(x,w(t))\ 2 



dx 



< 



JCl 



-2e 



dw , . 



dxdt+ {— + Y 



esssup / x|B(a;, w(t))\ 2 dx 
te[o,r] Jn 

+ [ X^~M^^o)\ 2 dx. 
Jn 2ei 



Choosing here e = 1/(2t) and = 1/2 we derive 

JessBup||xB(a!,ti;(*))||^n<T xD(z, ^) * „ , 
^ te[o,T] ^ ct / 2,f2x(o,T) 



+ ||xl>(a:,«7o)|llo, (4.13) 



Vr e [0, T]. We end up with estimate (|4.1ip . by combining inequality (I4.13[) with 
energy estimate (14.101) and Korn's inequality 



\(j>\\ o <£7 fc (n)||D(a:,^)|| 2 , n> 



(4.14) 



which is valid for all functions <fi e L 2 (S1) such that D(x, 0) € L 2 (f2) and vanishes 
in the trace sense on some open subset of <9f2 [HJ Chap. Ill, Sec. 3.2]. □ 

o 

Corollary 4.5. Let F, * e L 2 (Q), v o ,0 o G £ 2 (^), and w <E_W\ {tt). If Model A 
is solvable, then there exists exactly one generalized solution. 

Proof. Since Model A is linear, the uniqueness assertion amounts to the propo- 
sition that, if F , vq, Wq = and "f, 9q = 0, then there is only trivial solution. The 
latter proposition is obvious due to energy estimate (|4.10[) . □ 

5. Galerkin's Approximations and Existence of Solutions to Model A 

o 

Let {0,} C C§°(n) 3 and {tpi} C C§°{n) be total systems in W£ (fi) 3 and 

o 

W\ (fi), respectively. 

We construct Galerkin's approximations of the displacement vector and of the 
temperature distribution in the forms 



w n (x,t) = ^a/(t)0 z (x), 8 n {x.t) = y^bi(t)ipi(x), 



(5.1) 



i=i 



i=i 
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where unknown functions ai(t) and bi(t) (I = l,...,ri) are found from Galerkin's 
system 

it, d I a rp{x)Mx) ■ 4>j(x)dx 

1=1 al Jn 

n A ( +\ f 

= -V] 37 / x(x)(a u div x cj) l (x)-div x cj)Jx)+a^B(x 1 cj) l (x)) :B(x,cj)Jx)))dx 
i=i dt J " 

71 r. 

-V»iC) / [x(x)a p div x cf) l (x) ■ div^O) + (1 - x(a;))a 7) div x 0,(x) • div^^x) 
1=1 J n 

+ (l-xix jjaxDix,^)) : B(x, 0,(x))]dx 
+ y^fyft) / ae{x)xl)i(x)<i\v x <j>j(x)dx + / a F p(x)F n (x,t) ■ 4> J (x)dx, (5.2) 



j = 1, . . . , n, where F n is a given approximation of F such that 

F„eC°°(Q), F„ -» F in L 2 {Q) as n/ oo, (5.3) 

n Ah (+\ f Up 

V— 1 -t— / c p (a;)V>Kx)V>j(x)dx = - V& ; (i) / K(x)V x ipi(x) ■ W x ^j(x)dx 
i=i at Jn i=i 

+ / de o(x)(div x 4> l (x))tP j (x)dx + / y n (x,t)ipj{x)dx, (5.4) 

j = 1, . . . , n, where ^ is a given approximation of "J/ such that 

*„eC°°(Q), *„ -» * in L 2 (Q) as n / oo. (5.5) 

Since a T p(x) > 0, c p (x) > and the sets (cf> l7 . . . , cf> n ) and (tpi, . . . ,tp n ) are 
linearly independent in L 2 (fl) 3 and £ 2 (f2), respectively, then the matrices 

A„ = ( / a T p{x)4> l {x) ■ 4>Jx)dx) , B„ = ( / c p (x)^ ; (x)^j(x)cZx) 

are invertible, due to the classical theory of Hilbcrt spaces. 
Hence, setting 

ci(t) = -^7^. a n{t) = (ai(t), a n (t)), 
K{t) = {bi{t),...,b n {t)), c n (t) = (ci (*),..., c(t)), 

A " } = {f n ^ x ^ ai/dWx(l> ^ ' div x 0,(x) + a (U ©(a;,0 i (x)) : ©(a;, ^.(x)))^)™^, 



Al 2 ) 



(y [x( a; ) Q! p d i v x0;(a;) • div x 0.,(x) + (1 - x^a^divx^x) • div x 0.,(x) 

)n 

A i 3) = ( / ae(a;)^(a;)div a; 0,-(a;)da; > ) 
B« = ( / x(x)V x VK^) ' V x ^j(x)dx) n , 
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lffi(2) 



( / ae(x)(div x <pi(x))ipj(x)dx] 

K Jn '1,3=1 

F n {t) = (^J a F p{x)F n (x,t) ■ d) j (x)dx 

*n(*)=( / Vn(x,t)lpj(x)dx) 



3 = 1 

we see that system (|5.2j) - (|5.4j) is equivalent to the system of the first-order linear 
differential equations with constant coefficients in the normal form: 

~Tj, = ^ c n{i) ~ A n hS n ^ CL n (t) 



'&„(*) +A" 1 F„(t) ) (5.6) 

da n (t) 



= -B- 1 B( 1 )6„(t)+B- 1 (B( 2 ))*c„(i)+B- 1 *„W, (5.7) 



dt 



= c n (t). (5.8) 



System (|5.2|) - (|5.4|) or, equivalently, system (|5.6|) - (|5.8p . is supplemented with 
initial data 

ai(t)\ t =o = °°i d - / w (x) ■ 0j(a;)dx, (5.9) 
Jn 

6,(i)| t=0 = &? = f / 6 (xWi(x)dx, (5.10) 



q(*)| t =o - c° / v (x) ■ 4>iix)dx. (5.11) 
Jn 



On the strength of the classical theory of systems of first-order ordinary lin- 
ear differential equations, the Cauchy problem (|5.6|) - (I5.11|) has a unique infinitely 
smooth solution (a n (t), b n (t), c n (t)) for any neff. This amounts to the following. 



Proposition 5.1. Galerkin's system (|5.2p - (|5.4j) . supplemented with initial data 
(|5.9p - (|5.11|) . has a unique smooth solution (ai(t), . . . , a n (t) , b\(t) , . . . , b n (t)) on [0, T] 
for any n e N. 

The approximate distributions of pressures now can be found from the equations 

p n (x, t) = -x( a; )apdiv x u;„(a;, t), (5-12) 

q n (x,t) =p n (x,t) + — , (5.13) 

a p ot 

n n (x,t) = -(1 - x{x))a v div x w n (x,t). (5.14) 

Remark 5.2. Using the standard technics [H Sec. III. 3] we easily conclude that 
five approximate functions (w n , 9 n , p n , q n , n n ), which are obtained by virtue of 
Eqs. (|5.1[) - (|5.4[) and (|5.9[) - (|5.14[) . are a generalized solution of Model A, provided 
with the given approximate functions F n (x,t) and ^f n (x,t), and with the initial 
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data 



w n (x,t)\ t=0 = ^a^tb^x), 
dw n {x,t) 



1=1 



at 



|t=o = ^2^cj) l (x), 
1=1 

n 



(5.15a) 
(5.15b) 
(5.15c) 



i=i 



where a?, 6°, and c° are given by Eqs. (|5.9p - (|5.11[) . i.e., they are Fourier coefficients 
of initial data Wq, v , and Qq, and hence initial data (|5.15[) are the partial Fourier 
sums of w Q , v , and 9$. 

Due to this remark and energy estimate (|4.10p . the sequence (w n , 9 n , p n , q n , 7r n ) 
has a weak limiting point (w, 9, p, q, ir) as n f 00, and, due to linearity of Model 
A, the functions w, 6, p, q, and ir are a generalized solution of Model A, provided 
with initial data w , v , and 9 , which completes the proof of Theorem 13.41 

6. Additional Estimates for the Pressures 

Energy inequality (|4.10[) includes some estimates for the pressures p, q, and ir, 
due to Eqs. (|3.6d[) - (|3.6fj) . However, these estimates are not applicable for analysis 
of incompressible limiting regimes, because they are not uniform in a p and a n , as 
one or both of these coefficients grow infinitely. Hence, in order to do such analysis, 
it is necessary to obtain additional bounds on the pressures. 

We start with justification of the following technical result. 



Lemma 6.1. Let F, dF/dt, dm /dt £ L 2 (Q), and initial data be homogeneous, 
i.e., 

w = 0, v = 0, 6» = 0. 
Then the following bound is valid: 

-09 



(6.1) 



1 

— ess supa T 



te[o,- 



d 2 w 2 

V?— (*) 



dt 2 



2.9. 



— ess sup a v 
2 te\o,r] 



-ess supct\ 
2 te[o,r] 



1 

— ess sup 

2 te[o,r] 

dw 2 
(1 - x)div x — {t) 



2 

2.9. 



2 si 



1 

- — ess sup a p 



te[o,T] 



_ dw 2 
Xdiv x —(t) 



, _ / dw , . 
(l-Xp(x,—(t) 

/ d 2 w\ 2 



2.9 



Xdiv x 



d 2 w 



dt 2 



2.0 



2,f2x(0,r) 



< 







r=dF 2 

2 r=dF 2 

VP 



1 



2a 



dt 



2,fix(0,r) 

H 

2,!)x(0,t) 

2,six(o,t) ^ 2 min{cp/, c ps } 



5 2 

V* ir yx9 

dt 



2 min{ c p /, c ps } 
1 



2,flx(0,r) 

dm 2 

2,nx(0,r) 



dm 


2 1 


dt 


2,Ox(0,i). 



1 



a r mm{ Cp/,c ps | 



*lt= 



t=o||2,n 



= CSW, Vre[0,T]. 

(6.2) 
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a T p- 



t=0 



(6.3) 



Proof. We take advantage of the fact that the generalized solution is unique 
and may be constructed, using Galerkin's method. 

First of all, we notice that the values of all the constants a?, &?, and c?, defined 
in Eqs. (|5.9p - (|5.11[) . are equal to zero. Passing to the limit in the right hand sides 
of Eqs. (or, equivalently, of Eqs. (pT6 |) -([0 |) ). we conclude that 

d 2 w n _ _ d9 n 

= a F pF n \ t=0 , Cp—- 
t=o at 

Next, differentiating Galerkin's system (|5.2 p ~(|5.4 p and equations (|5. 12[) (15. 14[) 
with respect to t, on the strength of Remark 15.21 and Eqs. (|6.3|) we conclude that 
the derivatives 

f dw n dO n dp n dq n dir n 

KdT' ~~dt> ~dt' ~dt : ~dt 
are the generalized solutions of Model A, provided with the dimensionless approx- 
imate density of distributed mass forces dF n /dt and volumetric density of heat 
application d^ n /dt, and with the initial data (|6.3|) and (|5.15b[) . Now passing 
to the limit as n /* oo in energy estimate (I4.10[) . we immediately derive bound 
(|cT2l . □ 



Next, we introduce the normalized pressures p, q, and tt as follows: 

v( x ) 

p{x,t) = -x(x)a p div x w(x, t) + 



meas ft f j n 
Tr(x,t) = -(1 - x(x))a v div x w(x,t) 
, 1 - X(x) 



q(x,t) = p(x,t) + 



meas Q s Jq 

a v dp(x,t) 



x(x)a p div x w(x,t)dx, (6.4) 
(6.5) 

(1 — x(cc))a^div a; w(a;, t)dx, (6.6) 



dt 



(6.7) 



Remark 6.2. By the straightforward calculation, using Green's formula, we con- 
clude that integral equalities (|3.8[) and f|4. 3|) with Eq. (|3.6bp being inserted into 
them, are valid with q and 7f on the places of q and 7r, respectively, and that 

pdx — I frdx = / qdx = 0. (6-8) 
n Jn Jn 

This means that all the results previously obtained in this article, in particular, 
Theorem l3.4l and energy estimate (|4.10p . remain true for the modified model, which 
appears if we substitute equations (|3.6dp - (|3.6fp in the statement of Model A by 
equations (j6T4|H(f6~7|) . 



For the normalized pressures p, q, and n we prove the following. 

Theorem 6.3. Let F, dF /dt, vp, /dt E L 2 {Q) and initial data be homogeneous, 
i.e., satisfy Eqs. (|6.ip . Then the normalized pressures p, q, and tt satisfy the 
estimates 



19111.0+11*11!,, 



(" 



l2,Q 



1*1 



2,Q 



dF 


2 


a* 


dt 


+ 

2,Q 


dt 



in=olli,n + ll*|t=o||l n 



def 



(6.9) 
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\\p\\l Q <c; nc + c en (T), 

where C inc = C inc (T, fl, a T , a F , a M , a\, agj, a Ss , pf, p s , c p f, c ps ). 



(6.10) 



Remark 6.4. We emphasize that the constants C en (T) and Ci nc do not depend 
on a p and a v , which implies that the obtained bounds are applicable for studying 
asymptotic as a p and a v grow infinitely. Moreover, notice that the constants C en (T) 
and Ci nc also do not depend on geometry of flf and fl s . This fact seems to be very 
useful in view of possible studies of homogenization topics in periodic structures 
(like in [I US]) for Model A. 

PROOF. On the strength of Lemma T6. 11 the bound 



VP 



d 2 w 



< 



dt 2 

T r 



2 1 

2,£1x(0,r) 2 

OF 2 



-de 

Cp dt 



+ 



2a T 

o 

e T - 



VP- 



dt 



2,f!x(0,T) 
1 



2,Ox(0,t) 



r=dF 2 



2min{c p /, c ps } 
1 



a* 


2 1 


dt 


2,nx(o,t). 



'dt 



2,ox(o,r) 2 min{cp/, c ps } 

2 1 



< 



VP 



VpF\t=o 
OF 2 



2,n 



+ 



2,Ox(0,t) 
1 1 2 



dt 



07 



+ ^\\VpF\t=o 



dt 

i2 

l2,n 



min{c p /, c ps } 
1 



t=0 



2,f!x(0,r) 
1 



mm 



in{c p /, Cp S } 



mm{c p f,Cps} 



t=0\ 



2S> 



2, SI 
d^J 2 

2,S1x(0,t) 

Vtg[0,T] (6.11) 



0f 



is valid. 

Justification of estimates (|6.10p - (|6.9p is based on use of this bound, on a special 
choice of test functions in integral equality (|4.3p . and on application of energy 
inequality (|4.10[) and Remark 16.21 

We integrate the first term in Eq. f|4.3|) by parts, substitute Eqs. l|3.6b[) and 
(|3.6el) . and replace p and n by p and tt, respectively, which is legal due to Remark 
Thus we get 



f ( _d 2 w _ ( 

Jnxio.r)^ dt 2 V 



dw 



IQx(0,t) 

(1 - x) a \^( x > w 



dt J 



D(x,ip) - a e 9div x ip - a F pF ■ ipjdxdt = 0, Vr e [0,T]. 

(6.12) 



def 

Here we denote z = q + tt for briefness. 

Now our aim is to choose a test function ip in this equality such that div x ip = z 
and all integrals make perfect sense and admit estimates in terms of ||5|| 2 ,r2x(o.-r) 
independently of a p and a n . Evidently, if this choice of cp is possible then it leads 
directly to a bound on ||5||2,ox(o,t)j uniform in a p and a n . We succeed to pick up 
such a test function as follows. Introduce successively ipi — (fi(x, t) as the solution 
of the Dirichlct problem for Poisson's equation on Q, for t <G [0,T]: 



A x ^i = z, tpi\en = 



(6.13) 
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and (f 2 — <p 2 (x,t) such that 

div x (p 2 = 0, ¥> 2 |an = V x y>i|an. 



(6.14) 



Note that due to the classical theory of elliptic equations, a solution of problem 
(|6.13|) exists, is uniquely defined by z, and admits the bound [H Sec. II. 7] 



Ml . < Ci(fi)||z|| 2 . n Vte[0,T], 
VKf (n) 



(6.15) 



and that, on the strength of the construction in Chap. I, Sec. 2.1] and property 
(|6.8|) . a vector-function cp 2 satisfying Eqs. (|6. 14|) may be found such that 

(6.16) 



\\f2Wwm < C , (fi)||V a! vi||2 1 en Vte [0,T], 
which immediately implies the bound 

Wf2\\w}(n) < C 2 (n)\\<pi\\ o ^ Vte [0,T], 



(6.17) 



thanks to the well-known trace theorem 4, Sec. 1.6]. In inequalities (|6.15p - (|6.17p 
the constants Co(f2), Ci(O), and 6*2(^1) depend merely on regularity of dfl and 
make sense, for example, for C 2 -piecewise smooth surfaces dil. 

The sum tp = V x <fi — (f 2 is a valid test function for integral equality (|6.12p . 
Inserting it and rearranging terms, we obtain the equality 



z 2 dxdt 



fix(0,r) 



/ _d 2 w / Ow\ 



<9lO s 



Qx(0, 



dw 

~dtJ 



l(x, (p 2 ) + (1 - x)a A D(a;, to) : D(x, V^i) 

- (1 - x)aAB(a;, to) : y> 2 ) - a^z - aF/5-F • (V x <^i - ¥> 2 )) da;c ^ ( 6 - 18 ) 

Vr e [0,T], due to Eqs. (j6~T3| and (|6~14| . 

Applying Holder's and the Cauchy-Schwartz inequalities, bounds (|6.11|) and 
(|6.15| - (|6.17p . and energy inequality (|4.10p to the right hand side of equality (|6.18p . 
we estimate 



l^lla,nx(o,T) - 



1/1 



1 



2\e! e 2 



VP 



d 2 w 



dt 2 



2,f2x(0,r) 



aF\\VP F \\ 2 2,nx(0,T) 



-(— — 

1 „ r- 



dt J 



2,fix(0,T) 



+ "a||(1 - Xp(x,w)\\t,nx(o,r) 



2e 5 

max{a M , a\}e 3 



a y ll2,f2x(O,T) 



t T + of) max{p/, p s }ei 



.._ , _ ,..9 max{a„, a^l-e^ , 

||B(a! > V a pi)||l inx(0jT) + g HP(^2)ll2,Ox(0,r) 



I|l2,f2x(0,r) 
2 



(a- + a F ) maxjpf , p s }e 2 



m&xia gf ,a gs }e 5 „_„ 
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< 



.Si £2 



r=dF 2 



1 



1 



2,f2x(o,r) min{cp/, c ps } 

2 



Of 



2,fix(0,- 



+■ ^\\VpF\t=o\\l a + r^ n{( ^ fjCps} \\*\t=°\\ln + aF\\VpF\\ln, { o,r) 



-(— — 

2U3 + £4 



1 



^llV^llU(°.r) + minfe/><v>} ll^ll a .nx(D,r) 



\VpF\\l,nx(o,t) + ^-7— — T 11*11 2,nx(o,t) 



1 max{g 8/ ,g fe } 
2e 5 min{c p /, c ps } 

a 2 



min{c p /, Cp S } 

^ll^ll^x(0,r) + min{Cp 1 /!Cps} 

ll*l' 2 



II Xj/ lll,£lx(0,-r) 



1 



^IIVp^lll,ax(o, t ) t min{Cp/)Cps }ll-ll2,nx(o,t) 
2 



e T - t dt 



(a T + a F ) max{/9/, Ps}(ei + C 2 (f2)e 2 ) 



+ max{ap, a x }(e 3 + C 2 (Q)£ 4 ) + max{a e /, a S s}£5 PH 2 , 



!)x(0,t) 



(6.19) 



Choosing 



£1 = 2/(hCl(VL){a T + a F )m&x{p f ,p s }), 
£ 2 = 2/(5C^Q)Ci(fl)(a T +« F )max{ P/ ,ps}), 
£ 3 = 2/(5C 2 (r>)max{a M ,a A }), 
£4 = 2/(5C7 2 (f>)C 2 (r!) max{a M , aA }), 
£5 = 2/(5C 2 (fi)max{a e /,aes}), 
we deduce the estimate 



?n 2 <a (1) 



"Fll /=r.||2 1 ||,T,I|2 



Iv^lli, 



nx (°^ + miii{c p/ ,c ps } 11 W ll2.nx(o,r) 



+ G 



(2) 



^<9F 2 



1 



-||^|t=o 



2,f2x(0,T) minjcp/, Cps} 
1 



9* 



0t 



2,!1x(0,t) 



^ - ndnfa,,, cps} "^'^"^ + a 'II^H2^(o, T ) 



(6.20) 



where 



CP = (5/4)(e r + l)Cf{n){l + C7 2 (f!)) max{ Q „ a A } 

+ (5/4)(e r + l)C 2 (f7)(max{a 0/ ,a e s}) 2 / min {cp/,Cps}, 

d 2) = (5/4)(e r + l)C 2 (ft)(l + C 2 (f!))(a r + a F )max{p f , Pa }. 
Since the supports of q and tt do not intersect, we have 



l^ll2,nx(o,T) 



IkT 



H2,Qx(0,r)' 



I2,Ox(0,t) 

Thus bound (|6.9p immediately follows from inequality (|6.20|) . Finally, estimate 
(|6.10| follows from bound ()6.9|) . energy estimate (|4.10p . and equation (|6.7jl . □ 
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7. Additional Estimates for the Deformation Tensor 

Also we are interested in investigating limiting regimes arising as ax grows in- 
finitely. In order to fulfill this study, it is necessary to establish additional bounds 
on solution of Model A, which should be in certain sense uniform in a\. 

In this line, the fundamental estimate immediately follows from energy inequality 
dHO]): 

12 ^ r < f_w„. w_ ^ rn ( 71 ) 



e SS su P \\(l-xMx,w(t))\\ 2 2M < C en (r)/a x , Vr e [0,T]. 
te[o,r] 



Remark 7.1. Assume that (1 — x)B(a;,w;o) = 0, which implies that C en (Y) is 
independent of a\. Then from bound ()7.1|) it obviously follows that 



||(l-x)D(ss,to(*))||a,n — >0 as a A / oo. 

Since the kernel of the operator <p i— > D(x, <p) is the set of absolutely rigid body 
motions, i.e., translations and rotations <fi(x) = x + lj x x (x q ,uj — const) [TH 
Chap. Ill, §2.1], we see that the infinite growth of ct\ leads to solidification limiting 
regimes. 

The following assertion provides one more useful bound in the case of potential 
(e.g., gravitational) mass forces under additional assumption on geometry of tt s : 

Assumption 7.2. Let the Lebesgue measure of <9f2 n d£l s be strictly positive and 
f2 s be connected. 

Theorem 7.3. Let Assumption \ 7.2\ hold. Suppose that initial data satisfy Eqs. 
(|6.ip . the derivative d"^ /dt belongs to L 2 (Q), and the given density of distributed 
mass forces has the form 

F(x,t) = \7 x $(x,t) (7.2) 

with some potential $ £ W$(Q) such that {d/di)V x § e L 2 {Q). 
Then the following bound holds true: 



-ess supa T a\ 
'■ te[a,T] 

+ a^ax 

1 a p 



-dw , 2 



2.S2 



+ OL v OL\ 



\div. 



dw 
~dt 



dw\ 
~dt) 



2,Q 



1 a v 

ess sup 

4 a\ t e [o,T] 



2.Q 

(1 - x)ce\div x w(t) 



2 

2.Q 



— — - ess sup 

4 Ct\ t g[o,T] 



Xa\div x w(t) 

* -( 1 + ^ + S)("*iW 

where C so i = C so i{T, fl, x, ag, olf, p, c p ) 



+ ~esssup\\(l-x)a\®(x,w(t))\\l n 

2 ^ 4 t g[o,T] 



dt x 



1*1 



2,Q 





2 ) 


dt 


2,Q/ 



(7.3) 



Remark 7.4. We emphasize that the constant C so i does not depend on a p , a v , 
and a\. Thus estimate (|7.3|) becomes uniform in ct\ /* oo, whenever it is assumed 
that ctri and a p also grow infinitely such that 



0(a x ), a p = 0(ct\) as a x / oo. 



(7.4) 
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Remark 7.5. On the strength of Eqs. (|3.6d|) - (13.6fj) and energy estimate (|4.10|) . 
we observe that 



a p 

—ess sup 

a A te[0,T] 



Xa\div x w(t) 



— L ess sup 

a A te[0,T] 



\q\\2,Q < \\P\\2,Q < \\p\\2,Q + a V 



(1 - x)a\div x w(t) 
dw 



Xdhv 



dt 



2,Q 
2 

2 : 
2.Q 



Ot\ / \ II 2 

— ess sup||p(i)|| 20 , 

a P ie[0,T] 

—ess sup ||7r(t)||| n , 

a r; te[0,T] 



< lbl|2,Q + y/a v C en (T). 



These formulas and estimate (|7.3p yield that the pressures p, q, and 7r stay 
bounded in L 2 (Q) as ot\, a v , and a p grow infinitely such that limiting relations 
(HH) hold. 



Proof of Theorem [TH Substitute Eqs. (|3~6bf . ([3T6d)) -((3T6f )> . (jfn|) into Eq. 
HOI) to get 



o Jn 



_dw d(p _ dw 

aTP ~dT ' ~dt ~ xa P dlVxW ' uWxP ~ Xa„aiv x — ■ div x <p 



dw \ 

Xa^Byx, —J : B(x, tp) - (1 - x)ar,div x w ■ div^y 

— (1 — x)aAlD ) (a;, u?) : D(.t, <p) + ae6div x tp — appF ■ ipjdxdt 
f duo 

= a T p—(x,T)-<p(x,T)dx, Vr€[0,T]. (7.5) 

On the strength of Lemma [6.11 tp = a\dw/dt is a valid test function for Eq. 
(|7.5p . Inserting it and Eq. (|7.2p into Eq. (|7.5[) . we arrive at the equality 



-a T a A 



VP^-( r ) 2n + ^ a Aa p ||xdiv x w(r)||^ n 



<9«; 2 

xdw.— 



2,f2x(0,r) 



(■ 



IF J 



2,f!x(0,T) 



+ ^«a||(1 - x)div !C w(T)||| )n + -||(1 - x)ax^{x,w(T))f 2>Q 

f T f ( dw dw \ 

= -J J(a F axpV g Q- — -a x a e 6div x — jdxdt, Vtg[0,T]. (7.6) 

Using integration by parts on the right-hand side of this equality, we get 



-a T a\ 



+ a u a\ 



2.n 



+ a v a.\ 
( dw\ 



dw 2 



2,f2x(0,r) 
1 



, ^a\a p \\xdiv x w(T)\\l n 

2,!)x(0,t) Z 



1 



+ sOtAOijIKl - x)div x ti;(T)||| ]n + ~||(1 - %)a A D(a;, w(r))|||j 



22 



A. M. MEIRMANOV, S. A. SAZHENKOV 



SUBMITTED TO EJDE 



= l aFOe\Pf§(x, T)div x w(x, r)dx — / / apaxpf-— div x wdxdt 
Jn Jo Jn ot 

- I a F a x (l ~ x{x)){p s - Pf)V x ${x,T) ■ w{x,r)dx 
Jn 

+ [ f a F a x (l - x)(Ps - Pf">(§i V ^) ' wdxdt 

+ / a\ae(x)9(x, r)div x w(x, r)dx 
Jn 

a x a s —div x wdxdt, VtG[0,T1. (7.7) 
a at 

Now, on the right-hand side apply the Cauchy-Schwartz and Korn's inequality 
(|4.14p (with £l s substituted for fi) along with a simple estimate 

ll0lll,nx(o,r) <resssup||0(t)||2 V <f> E L°°(0, T; L 2 (fl)) (7.8) 
te[o,r] 



to get 
1 



-ess supa T aA 
J te[o,r] 



i—dw , . 



2.n 



+ Ot v Ot\ 



Xdiv ; 



dw 



2,Ox(0,t) 



1 



dw\ 2 



dt ) 



+ -esssupa A a^||(l - x)<&Vxw(t)\\l n 

2,nx(0,r) 2 t£[o, T ] 

esssupa x a p \\xdiv x w(t)\\l >n + -esssup||(l - x)<xxB(x,w(t))\\l n 
1 te[o,r] 1 te[o,T] 



< 



2 2 

° 7 % SS sup||(l-x)$(t)|||o 



2e 



te[o,T] 



J -esssup||x$(i)||2,o 



2 2 



2e 2 



(i-x) 



At 



2ei 

2 2 



2,nx(0,r) 2£3 



*e[o,r] 



2,fix(0,r) 



~ Pff 

2e 4 



esssup||(l-x)V a; $(0|| 
te[o,T] 



2.<> 



^-^ 2 |(i-x)|v^ 2 



2e 5 



2,nx(o,- 



'^-esssup ||(1 - X)«(*)||l.n + ^esssup ||x<9(*)||2,n 



2^6 te[o,r] 



2e 8 



(i-x) 



00 



of 

2,nx(o,T) 2eg 



2£7 te[o,T] 



X 



dt 



2,f!x(0,r) 



-^ (eo + t£2 + £e + Te 8 )esssup ||(1 - x)div x w(t)\\l n 
2 te[o,r] 

+ -rr( £ i + r£ 3 + e 7 + re 9 )esssup ||xdiv x u;(t)||2 )n 

te[o,r] 

(e 4 + re 5 )ess sup || (1 - x)a1;D(x, «>(*))||2,n> 
te[o.r] 



a 
"2 

(c fe (a.)) 2 



Vre[0,T]. (7.9) 
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Choose 



So = £6 = , El = £7 = 7; ) £2 = &8 = -3 • 

a„ 1 1 



£3 £9 8ra A ' 64 4(C fe (^)) 2 ' £5 4r(C k (Sl s )) 2 ' 

and set t — T. Applying energy inequality (|4. 10[) and Lemma 16.11 in order to 
estimate the norms of 9 and d9/dt in (|7.9p . we arrive finally at inequality (|7.3p . □ 



8. Incompressibility Limits 

In this section we prove the theorem, which explains the limiting behavior of 
solutions of Model A, provided with homogeneous initial data, as the coefficients 
a p and grow infinitely: 

Theorem 8.1. Let F, dF/dt, /dt £ L 2 (Q), and initial data be homoge- 

neous, i.e., satisfy Eqs. (16. ip . Suppose that in Model A the positive constants pj, 
p s , a T , a v , a^, a\, agj, ao s , ctp, c p f, c ps , Hf, and x s are fixed and the coefficients 
a p and depend on a small parameter e > 0: Let 

lim at. = a° lim ai = a° (8.1) 

e\o p p s\o 1 n 

For a fixed e > by (w s , 9 E ,p 6 , q 6 , tt e ) denote the generalized solution of Model 
A such that the pressures p s , q £ , and n E are normalized, i.e., satisfy equations 
(HHD -dllTJ) instead of Eqs. j3Jd| -^M ■ 
Then the following assertions hold true. 

1) Lfa p S (0, oo) and — oo, then the sequence (w e , 9 e ,p e , q s , tt £ ) is convergent 
such that 

o 

weakly* in L°°(0, T; W\ (SI)), (8.2) 

O 

weakly in L 2 (0,T;W% (SI)), 

weakly* in L°°(0, T; L 2 (Sl)), (8.3) 
weakly in L 2 (Q), (8.4) 

and the five limiting functions (w,9,p,q,ir) are the generalized solution of Model 
Bl with a p — a p . ( Statement of Model Bl and the notion of its generalized solution 
are given immediately after formulation of the theorem.) 

2) Ifctp = oo and € (0, oo), then the sequence (w e , 9 e ,p e , q £ , ir e ) is convergent 
such that the limiting relations (|8.2| - (|8.4p hold true and, moreover, p — q and the 
four limiting functions (w,9,q,ir) are the generalized solution of Model B2 with 
a, q = a®. ( Statement of Model B2 and the notion of its generalized solution are 
given after statement of the definition of generalized solution of Model Bl.) 

3) If a p = a° = oo, then the sequence (w £ , 9 s , p e , q £ , ir e ) is convergent such that 
the limiting relations (I8.2p - (|8.4p hold true and, moreover, p — q and the triple 

def 

of limiting functions (w,9,Sj = q + n) is the generalized solution of Model B3. 
( Statement of Model B3 and the notion of its generalized solution are given after 
statement of the definition of generalized solution of Model B2.) 



w — > w, 



dw £ aw 



\o dt e\fl dt 



e\,0 



p — > p, q — >q, 7T — >tt 

eVO e\0 eVO 
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Statement of Model Bl. In the space-time cylinder Q it is necessary to find a 
displacement vector w, a temperature distribution 9, and distributions of pressures 
p, q, and ir, which satisfy equations (|3.6a|) - (|3.6c[) . (|6.4p . and (|6.7|) , incompressibility 
condition in the solid phase 

(1 - x)div x w = 0, (x,t)eQ, (8.5) 

homogeneous initial data (|6.1|) . and homogeneous boundary conditions (|3.6hj) . 

Definition 8.2. Five functions (w,9,p,q,ir) are called a generalized solution of 
Model Bl if they satisfy regularity conditions (13. 7|) . condition ir € L 2 (Q), initial 
data ()6.1|) . integral equalities ()3.8|) and (|3.9|) . and equations (I3.6bj) . (16.411 . (16. 7| . and 
([53)1 a.e. in Q. 

Statement of Model B2. In the space-time cylinder Q it is necessary to find a 
displacement vector w, a temperature distribution 8, and distributions of pressures 
q, and 7r, which satisfy equations (|3.6ap - (|3.6cp and (|6.6p . incompressibility condition 
in the liquid phase 

Xdiv x w = Q, {x,t)eQ, (8.6) 
homogeneous initial data (|6.ip . and homogeneous boundary conditions (|3.6h[) . 

Definition 8.3. Four functions (w, 9, q, tt) are called a generalized solution of Model 
B2 if they satisfy regularity conditions ()3.7j) . condition q E L 2 (Q), initial data (16. 
integral equalities (|3.8|) and (|3.9p . and equations (|3.6b[) . (|6.6p , and (|8.6p a.e. in Q. 

Statement of Model B3. In the space-time cylinder Q it is necessary to find a 
displacement vector w, a temperature distribution 0, and a distribution of pressure 
which satisfy equations (|3.6ap . (|3.6cp . and 

F=-(ft + a e 6)I + xa^x,—) + (l-x)a x ®(x,w), (x,t)eQ, (8.7) 
incompressibility condition in the both phases 

div x io = 0, (x,t)eQ, (8.8) 
homogeneous initial data (|6.ip . and homogeneous boundary conditions (|3.6hp . 

Definition 8.4. The triple (w^,^) is called a generalized solution of Model B3 
if it satisfies regularity conditions (|3.7p , condition S) € L 2 (Q), initial data (|6.ip . 
integral equalities (|3.8p and (|3.9p . and equations (|8.7p and (|8.8p a.e. in Q. 

Remark 8.5. Clearly Theorem 18.11 provides the existence results for generalized 
solutions to Models B1-B3, as a byproduct. At the same time, existence and unique- 
ness of generalized solutions of Models B1-B3 may be justified without consider- 
ing Model A firstly, but starting from Definitions 18.2118.41 introducing the proper 
Galerkin's approximations, and keeping track of the proof of Theorem 13.41 Thus, 
Models B1-B3 are well-posed. 

Remark 8.6. Eqs. (|8.5p . (|8.6p . and (|8.8p are indeed the conditions of incom- 
pressibility, since, on the strength of continuity equations (I2.25d[) and (|2.25hp . they 
imply that p = p s in fl s x (0, T), p = pf in Qf x (0,T), and p = p in Q, respec- 
tively. Consequently, Models Bl, B2, and B3 are respective solid, liquid, and total 
incompressibility limits of Model A. 
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Proof of Theorem 18.11 We verify assertion 3 of the theorem only. Justifica- 
tion of assertions 1 and 2 is quite similar. 

If a £ a £ — > oo then there exist a subsequence (w £ , 9 e ,p e , q £ , ir £ ) and a limiting 

set of functions (w, 9,p, q, it) such that limiting relations (|8.2j) - (|8.4j) hold true, due 
to energy estimate (I4.10j) . Lemma [6.11 Corollary 14.41 Theorem 16.31 and homoge- 
neous initial conditions (|6.ip . Energy estimate (|4.10|) yields that 

es SS up\\xdW x w e \\ 2 2n < 2(C en (T)/a £ ) — ► 0, 
te[o,T] £ N.° 

e SSS up\\{\ - x)&v x w s \\l n < 2(C en (T)/a e ) — ► 0, 
te[o,T] 

since initial data are homogeneous and therefore C en {T) does not depend on e. 
Hence Eq. (|8.8jl holds true for the limiting function w. 

Next, substitute Eq. (|6.4[) into Eq. (|6.7p and pass to the limit as £ \ 0, using 
limiting relations (|8.2p and (|8.4|) , to get 



q = p— x a u{d/dt)div x w. 

On the strength of Eq. (|8.8p this equality yields that q — p on Q. 

Now, on the strength of limiting relations (|8.2p - (|8.4p . passing to the limit in Eqs. 
(|3.8p and (|3.9p as e \ 0, we conclude that the triple (w, 9, q + it) is a generalized 
solution of Model B3. 

In order to complete the proof of assertion 3, it remains to notice that due to Re- 
mark 18.51 the solution of Model B3 is unique. Hence the sequence (w £ , 9 e ,p e , q £ , ir £ ) 
has exactly one (weak) limiting point and therefore converges entirely, and there is 
no need to shift to a subsequence. □ 

9. Solidification Limits 

In this section we observe limiting behavior of solutions of Model A, as the 
coefficient ot\ grows infinitely. 
First, we prove the following. 

Theorem 9.1. Let F, dF/dt, d^/dt £ L 2 (Q), and initial data be homoge- 
neous, i.e., satisfy Eqs. (|6.ip . Let Assumption \7.S\ hold. 

Suppose that in Model A the positive constants pf, p s , a T , a p , a v , a„, a M; agf, 
ag s , ap, c p f , c ps , Xf, and x s are fixed and the coefficient ct\ depends on a small 
parameter e > 0, ct\ — a\. Let 

lim a £ x = oo. (9-1) 

By (w £ ,9 £ ,p £ ,q £ ,tt £ ) denote the generalized solution of Model A corresponding to 
a fixed e > . 

Then the sequence {w £ ,9 £ ,p £ ,q £ ,tt £ ) is convergent such that 

w £ ^ w > weakly* in L°°(0,T;Wl (0)), (9.2) 

e\o at e\o at 

(1 - x)^xW £ — > 0, tt £ — > strongly in L°°(0, T; L 2 {Q)), (9.3) 

e\0 e\0 



e\0 



weakly in L 2 (0,T;W^ (fi)), 
weakly* in L°°(0, T\ L 2 (ty), (9.4) 
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p £ — > p, q £ — ► q weakly in L 2 (Q), (9-5) 

and the four limiting functions (w,9,p,q) are the generalized solution of Model CI. 

Statement of Model CI. In the space-time cylinder Q it is necessary to find a 
displacement vector w, a temperature distribution 9, and distributions of pressures 
p and q, which satisfy the equations 

d 2 w 

a T p f —2-=div x F f +a F pfF, mQ f x(0,T), (9.6a) 

F f = -(q + ao f 6)I + a^(x,^^, mQ f x{0,T), (9.6b) 

(1 - x)w = 0, in Q, (9.6c) 

09 Q 
c p — = div x (xV x 9) - xaef-^_div x w + in Q, (9.6d) 

p + xckpdiv x u; = 0, in Q, (9.6e) 

1 = P+— % in «> ( 9 - 6f ) 
a p at 

homogeneous initial data (|6.1[) . and homogeneous boundary conditions 

w = 0, for x G dCl f , t > 0; = 0, for x G dCl, t > 0. (9.6g) 

Definition 9.2. Four functions (w, 9, p, q) are called a generalized solution of Model 
CI if they satisfy the regularity conditions 

w, —,S7 x w,9,X7 x 9eL 2 (Q), (9.7) 
ot 

equations (|9~6bj) . ([9~6c]) . (f9T6ej) . and |9^6f|l a.e. in Q, and the integral equalities 

Jo J ( a rPf^f-^-- p f-^^ + a FPf F - ( p)dxdt^O (9.8) 
for all smooth y> = tp(x,t) such that = ^|f=r = and 

Jq ( dp9 l& ~ * Vx6 ' + Xae f (diy x w)^- + *^») dsu* = (9.9) 
for all smooth ip = ip(x,t) such that ip\dQ = ip\t=T = 0. 

Remark 9.3. Like in Remark 18.51 we notice that Theorem 19.11 provides the ex- 
istence result for Model CI as a byproduct, and that existence and uniqueness of 
generalized solutions of Model CI may be justified independently of Theorem 19.11 
by considerations similar to the proof of Theorem [331 



Proof of Theorem 19.11 Firstly we have that limiting relations (|9.2j) . (19. 4|) . 
and (|9.5|) hold true for some subsequence (w £ ,9 £ ,p £ ,q £ ) due to energy estimate 
P~Tu| and Lemma HU 

Secondly, bound (|7.ip immediately implies that the limiting displacement vector 
w satisfies the equality (1 — x)B(x, w) — in Q. On the strength of Remark 17.11 
and Assumption 17.21 this equality yields that Eq. (|9.6c| holds. From bound (|7.ip . 
Assumption 17.21 homogeneity of initial data (see Eq. I|6.1[) ). and Korn's inequality 
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(|4.14| (with fl s substituted for fi) it follows that 
ess sup (||(1 - x)w £ (t)\\ln + 11(1 - x)V^ £ (i)||i, n ) 

t£[0,r] 

< (Cfc^^esssuplKl-x)©^^^))!!^-^ ' 
te[o,r] £ \° 

which proves limiting relations (|9.3p . (By Cfc(f2 s ) the constant in Korn's inequality 
is denoted.) 

Next, on the strength of limiting relations (|9.2j) - (|9.5j) . passing to the limit in 
Eqs. (|3.8[) . (|3.9[) . (|3.6d|) . and (|3.6ep as e \ 0, we conclude that the four functions 
(w, 9,p, q) are a generalized solution of Model CI. 

In order to complete the proof, it remains to notice that due to Remark 19.31 the 
solution of Model CI is unique. Hence the sequence {w £ , 9 e ,p £ , q £ ) has exactly one 
(weak) limiting point and therefore converges entirely and there is no need to shift 
to a subsequence. □ 



Finishing this article, we derive a rather peculiar asymptotic, which is observed 
thanks to additional bound (|7.3p . Namely, we prove the following. 

Theorem 9.4. Let Assumption \7. 2\ hold, initial data be homogeneous, i.e., satisfy 
Eqs. (|6.ip . the given density of distributed mass forces have the form 

F(x,*) = V x *(x,*) (9.10) 

with some potential $ 6 W^{Q) such that (<9/<9i)V ;r < l > £ L 2 (Q), and the given 
volumetric density of exterior heat application be such that ^,d^/dt G L 2 {Q). 

Suppose that in Model A the positive constants pf, p s , a T , a v , a^, agf, ag s , 
a F> c p f, c ps , >r/ , and k s are fixed and the coefficients a\, a p , and a r; depend on 
a small parameter e > Q: a\ — a e x , a p = a e p , and a n — a^. Let limiting relations 
(|7.4p hold, or, equivalently, 

lima?=oo, lim ^ = a° e (0, oo), lim % = a° € (0, oo). (9.11) 

By {w e ,9 e ,p e ,q e denote the generalized solution of Model A corresponding to 
a fixed e > . 

Then the sequence (w E , 9 E ,p e , q £ , tt e ) is convergent and there exist functions 
(it, 9,p) such that 

w e — > 0, — > strongly in L°°(0, T; L 2 {tt)), (9.12) 

eNO at s\0 



e\0 



weakly in L 2 {0,T;W% (f2)), 

weakly* m L°°(0, T; L 2 (ty), (9.13) 

p £ —^p, q £ —^p, (9.14) 

tt £ -(l-x)a° v dw x u, (9.15) 

(1-X)a e x w e — >u weakly* in i°°(0, T; i 2 (fi)), (9.16) 

(1 -xp{x,a\w e ) — > (l-x)B>(x,u) weakly* in ^(O, T; L 2 (0)), (9.17) 

e\0 
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and the triple (u,0,p) is the generalized solution of Model C2, formulated below. 
Statement of Model C2. Find successively the following: 

( C2a) a temperature distribution 6, solving the initial-boundary value problem for 
the heat equation 

d6 

c p — = div x (xV x 6) + in Q, (9.18) 

0|t=o = 0, on SI, 9\ an = 0, for t > 0; (9.19) 

( C2b ) a hydrostatic pressure p in the liquid phase from the equilibrium equation 

p = -a sf 6 + a F p f $, in Q f x (0, T); (9.20) 

( C2c) an up-scaled displacement vector field u, solving the mixed problem for the 
3D system of stationary wave equations 

dw x B(x, u) + a°W x div x u = a Ss V x 6 - a F p s W x ^, in D. s , (9.21) 

B(x, u) ■ n + aP(div x u)n = {a es 9 - a F p f $)n, on T, (9.22) 

u = 0, on (dCl s \T). (9.23) 

(Here n is the outward unit normal to dQ s . Variable t appears in the system as 

def 

a parameter. Also, recall that T = dfl s n dilf and that the Lebesgue measure of 
90 s \ r is positive due to Assumption ^. 21 ) 

Definition 9.5. The triple (u,9,p) is called a generalized solution of Model C2 if 
it satisfies the regularity conditions 

9, V x 6,p, {1 - x)VxU £ L°°(Q,T;L 2 {n)), (9.24) 

and the integral equalities 



J ((XP- (1 - x)a°div x u + a 9)div x (p 



- (1 - x)D(ar, u) : D(x, cp) + a F pS7 x <$> ■ cp) dxdt = (9.25) 
for all smooth cp = cp(x,t) such that cp\gn = v\t=T = and 

/ ( dp9 l^ ~ * Vx6 ' Vxi> + yi>)dxdt = (9.26) 
for all smooth -ip = ip(x,t) such that ip\an = ip\t=T = 0. 



Remark 9.6. By the standard arguments it is easy to verify that integral equalities 
(|9.25p and (19.26)) are equivalent to equations (|9.18|) (|9.23|) in the distributions sense. 
Due to the well-known theory of second-order elliptic and parabolic equations [6l [7] 
there exists exactly one generalized solution of Model C2, whenever <!> and \& satisfy 
conditions of Theorem l9.4l and whenever we prescribe some certain values to (1— x)p 
and x u - (For example, we may simply set (1 — \)P — X u = 0.) 

PROOF of Theorem 19.41 On the strength of energy estimate (|4.10p . bound 
(17. 3|) . limiting relations (|9.1ip . and Remark l7.5i we conclude that limiting relations 
(|9.13| - ()9.17p hold true, at least for some subsequence, and that limiting relation 
(|9.12|) and relations 

XQ.div,— >0, xa^[x,—)^0 strongly in L 2 (Q) (9.27) 

OX e\,0 V Ot J e\,0 
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are fulfilled for the entire sequence e \ 0. 

The limiting functions u, p, and 9 satisfy regularity conditions (I9.24| due to 
relations ()9.13[1 — (|9.17[) and Korn's inequality 

K*)llwJ(n.) < C fc (fi.)||D(x,u(t))|| 2 ,n., Vt e [0,T]. 

Substituting Eqs. ([3lE]) -(|3lif ]) . and $LTU\i into integral equality (O, and 

passing to the limit as e \ 0, we arrive at integral equality (|9 . 25[) due to relations 
(I9.12p - (|9.17p and (|9.27[) . Analogously we justify that 8 satisfies integral equality 
(|9.26l) . Finally, it remains to notice that, since solution of Model C2 is unique, the 
sequence ((1 — x)a|u> £ , 9 £ , p £ , q e , ir e ) has exactly one (weak) limiting point and 
therefore converges entirely. □ 
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